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^h! Abstract 

Q-i, We study final-state interactions in nonleptonic weak decays in 

| statistical models by averaging over ensembles of strong interaction S- 

matrices. The models range from one with completely random strong 
• rH , interactions, which gives extensive mixing between physical states, 

r> | to models with feeble final-state interactions, characterized by small 

phase shifts. We compute expectation values for weak decay rates 
and the fluctuations around these means. The coherence between 
interfering amplitudes is gradually destroyed as the strength of the 
final-state interactions is increased. Provided there are at least two 
weak phases, there is direct CP violation as a result of the strong 
phases from final state interactions. Data indicate that for D meson 
decays, models with extensive mixing are appropriate, while for B 
meson decays, models with small final state interactions are required. 



*This work was supported in part by the Director, Office of Energy Research, Office of 
High Energy and Nuclear Physics, Division of High Energy Physics of the U.S. Depart- 
ment of Energy under Contract DE-AC03-76SF00098 and in part by the National Science 
Foundation under Grant PHY-95- 14797. 



Disclaimer 



This document was prepared as an account of work sponsored by the United States 
Government. While this document is believed to contain correct information, neither the 
United States Government nor any agency thereof, nor The Regents of the University 
of California, nor any of their employees, makes any warranty, express or implied, or 
assumes any legal liability or responsibility for the accuracy, completeness, or usefulness 
of any information, apparatus, product, or process disclosed, or represents that its use 
would not infringe privately owned rights. Reference herein to any specific commercial 
products process, or service by its trade name, trademark, manufacturer, or otherwise, 
does not necessarily constitute or imply its endorsement, recommendation, or favoring by 
the United States Government or any agency thereof, or The Regents of the University of 
California. The views and opinions of authors expressed herein do not necessarily state or 
reflect those of the United States Government or any agency thereof, or The Regents of 
the University of California. 



Lawrence Berkeley Laboratory is an equal opportunity employer. 



ii 



1 Introduction 



The study of CP violation is, to a large extent, the study of phases. Four- 
fermion interactions introduce CP violation if their coefficients cannot be 
made real through redefinition of the particle fields. Kobayashi and Maskawa 
[[|] recognized that, for three generations of quarks, mixing would in general 
leave an irremovable phase, guaranteeing the presence of CP violation, with- 
out the introduction of any new particles (aside from the c, b, and t quarks, 
which were then yet to be discovered!). It is not only through weak inter- 
actions, however, that phases are introduced into the amplitudes of weak 
decays. When the final states include hadrons, and in nonleptonic decays 
in particular, there are strong final state interactions that produce phases. 
Strong phases arise in accordance with unitarity, which requires imaginary 
parts for amplitudes of processes that have on-shell physical intermediate 
states. While the phase of a bare weak amplitudes changes sign when the 
process is taken to its CP conjugate, the final-state interaction phase is un- 
changed, reflecting the CP invariance of strong interactions. 

Final-state interaction phases can be either a blessing or a curse. Some 
potential CP violations are observable only if there are final state interactions 
present, with the measurable quantities being proportional, typically, to the 
sine of the difference of two strong-interaction phases. In other situations, 
the final-state interactions obscure the phases of the KM matrix that one 
wishes to measure. Both circumstances demand understanding of final-state 
interactions. 

Unitarity and time-reversal invariance (of the strong interactions) deter- 
mine final-state interaction phases. In the familiar case of K° — > tttt there are 
only two open channels, 1 = and 1 = 2. The weak decays pick up phases 
exp(i5j), where 5j is the strong-interaction phase shift in the appropriate 
isospin channel. Unfortunately, this simplicity is lost once there are many 
channels open to the strong interactions, as is the case for D and B decays. 
Formally, it is easy to describe the effect. If there are N open channels, we 
diagonalize the S-matrix to find N eigenchannels, each of which has an asso- 
ciated eigenphase-shift. We then analyze the weak decays in this eigenbasis 
and attach the phase exp(i5 a ) to the eigenchannel \a >. 

The difficulty here is two-fold. We do not know a priori the eigenphase nor 
do we know the composition of the eigenchannels. These problems involve 
the full complexity of non-perturbative physics in a particularly awkward 
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way. In this paper, our approach is statistical, based on the assumption that 
when the number of open channels, N, is large, it may be possible to study 
the general behavior of the decay amplitudes, even when it is not possible to 
make a precise prediction for any particular one of them. 

How big are final-state interactions in the decays of heavy mesons? At 
one extreme, it might be argued that they should be small. An appeal can 
be made to asymptotic freedom, as has been argued, for example, by Bjorken 
0. An explicit implementation of this philosophy was already given by Ban- 
der, Silverman, and Soni in 1979 ||, who calculated a one-loop diagram (a 
"penguin") that gave an imaginary part to a decay amplitude, using pertur- 
bative QCD. At the other extreme, one can cite the highly inelastic nature of 
strong collisions at high energies and the success of models with absorption 
in low partial waves ||]. 

A crude but simple measure of the extent of final-state interactions is to 
look at decays to final states that involve interference of two distinct isospin 
channels. Consider, for example, D° — ► K~tt + , K°tt°, whose final states are 
superpositions of I = 1/2 and / = 3/2. The "factorization" model argues 
that the charged current in weak decays should produce a charged 7r more 
frequently than neutral ones, by a factor 18 (even before QCD corrections, 
which further increase the dominance), requiring coherence between the I = 
1/2 and 1 = 3/2 amplitudes. In fact, the charged ir decay has a branching 
ratio just a factor of two higher than the neutral 7r decay. On the other hand, 
the decay B° -> D~ ~p + is more than an order of magnitude more probable 
than B° — > D p°, consistent with factorization. If final-state interactions 
in the 1 = 1/2 and 1 = 3/2 channels were substantial the relative phases 
between the I = 1/2 and 1 = 3/2 would have been randomized, making a 
large discrepancy between the branching ratios improbable. These examples 
suggest that at the higher energies of B decays, final-state interactions are 
less important than at the lower energies of D decays. 

Below we discuss a model that enables us to interpolate between the 
extremes of complete final-state interaction mixing and the complete absence 
of such mixing. 

Our point of departure is a statistical model of strong interactions based 
on random S-matrices as introduced (initially for nuclear physics) by Dyson 
|J. Dyson showed that each S-matrix can be expressed in terms of an N- 
dimensional orthogonal transformation (rotation) and a set of eigenphase 
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shifts. The weight attached to each S-matrix is a product of the weight 
associated with the rotation group and a weight associated with the particular 
set of eigenphases. This very elegant development is awkward to use, so we 
consider some simpler alternatives, as well. The random S-matrices lead to 
nearly complete mixing of all the communicating channels. This may be 
appropriate to the D decays, but is not for the B decays. We develop an 
alternative in which the phase shifts are restricted to a smaller range. 

2 Summary of Results 

The statistical approach of averaging over possible final-state interactions is 
conceptually straightforward, but it does involve some tedious calculations, 
as described in subsequent sections and as displayed more fully in some of the 
appendices. The results of these calculations, however, are simple to grasp 
and the essential points are displayed in this section. 

We imagine a weak decay amplitude X b to a physical state b in the absence 
of strong interactions corrections. The decay rate into b is then, up to a fixed 
constant and phase space factors, 

\M b \ 2 = \X 2 b \. (2.1) 

In the presence of final-state interactions, a decay that "initially" pro- 
duces a state a can rescatter into the state b. The extent of this rescat- 
tering can be modeled by introducing a distribution of phase shifts that 
are symmetrical about zero (so < sin 5 >= 0). The strength of the phase 
shifts is characterized by < cos 5 > and < cos 25 >. In the extreme of full 
final-state interactions, all values of the phase shift are equally likely and 
< cos 5 >=< cos 25 >= 0. In this instance, we find 

< \M b \ 2 >= ^(2|X b 2 | + NW\), (2-2) 

where 

W\ = ^E\Xa\ (2-3) 

is the mean of the absolute squares of all the "bare" weak decay amplitudes. 
For large N, most of the memory of the initial amplitude to b is lost, and 
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instead it is only the average of the square of the amplitude that counts. 
More generally, if smaller values of the phase shifts are favored, we find 

< i*ur > = N 



N + 2 
N 

N + 2 



-\X b \ 2 + \X 2 \ + [\X b \ 2 - \X 2 \) <cos5>' 



^< cos5 > 2 +j^)\X b \ 2 + (1- < cos5 > 2 )\X 2 \ 



(2.4) 



We can see how this interpolates between the situation with no final-state 
interactions, where only the amplitude to b matters, and the fully mixed case. 
We can study not only the expected means of the squares of the amplitudes, 
but also the fluctuations about those means. For large final state interactions, 
the fluctuations are comparable to the mean squared amplitudes themselves. 
For small final-state interactions the fluctuations of course are small. 

The coherence or incoherence of amplitudes for final states with different 
isospin is an important feature of weak decays and reveals much about the 
decay dynamics. In the absence of final-state interactions, the actual decay 
amplitudes are simply the bare decay amplitudes and we can write 

< \M b + M v \ 2 >= \X b + X b ,\ 2 . (2.5) 

Here the primed and unprimed amplitudes correspond to two states that 
cannot mix through strong interactions, because of differences in isospin or 
some other good quantum number. Introducing final-interactions mixes the 
amplitudes with other amplitudes, but with the primed and unprimed groups 
remaining separate. If there are N unprimed states and N' primed states, 
we find 

< \M b + My\ 2 > 

= I < cos 5 > X b + < cos 5' > X b ,\ 2 
W(l-<cf>') (1, x! , +W 



N + 2 \N 



N'(l- < cos 5' > 2 ) ( 2 , ... 2 



+ N> + 2 [WW + |X ' 2| J 

(2.6) 

If the final-state interactions are not large, we can expand in the phase 
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shifts and obtain 

<\M b + M b '\ 2 > — > | < cos5 > X b + < cos5' > X b ,\ 2 

+ <5 2 > ]X T \+ < 5' 2 > \X^\ + 0{l/N) (2.7) 

We see that the rate receives contributions from a purely coherent piece and 
from a purely incoherent piece, with the balance determined by the extent 
of the final-state interactions. There is complete coherence in their absence 
and complete incoherence when the phase shifts are completely random. 

Final-state interactions are essential to direct CP violation. We find, to 
leading order in 1/N, that the expected size of such CP violation is given by 



< (\Ml\ - \M 2 b \f >= ^ E |X a | 2 |X c | 2 sin 2 (0 a - C )(1- < cos 26 > 2 ). 

a,c 

(2.8) 

We see the expected structure. CP violation requires at least two different 
weak phases, a ,0 c , together with some final state interactions. If the 
phase shifts are small, so are the direct violations of CP. 

A comparison of some of the predictions with extant data on D and B 
decays is given in Section [TO, to which the anxious reader can skip. 



3 Unitarity and Final-State Interactions 

We are interested in decays of spinless particles so we restrict our considera- 
tion to states of zero angular momentum. It is conceptually simpler to think 
of the final states as two-body so that they are discrete, say N of them. The 
strong-interaction S-matrix in a fixed partial wave is thus N x N and in fact 
it can be taken to be symmetric, as a consequence of T invariance. It satisfies 
the unitarity relation SS^ = I. 

The symmetric S-matrix can be diagonalized by an orthogonal transfor- 
mation: 

S diag = OSO T , (3.1) 

where Sdiag has diagonal elements e 2lSa ; a = 1, ...N. The real N x N matrices 
O satisfy 
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o- 1 = o T , 

E O aa O a b = S a b, 
a 

E aa 0p a = 8 a p, 

a 

(3.2) 

where the Greek indices label states in the eigenbasis, while Latin indices 
label physical states. 

As shown in Appendix [A], the weak decay of a particle into the eigenstate 
a has an amplitude of the form 

M a = e i5 <*e i4 "*W a , (3.3) 
while the the CP conjugate decay has the amplitude 

Ma = e i5a e- i(t "*W aj (3.4) 

where W a is real, that is, the weak phase, 0, changes sign, while the strong 
phase, 5, does not. The decay amplitude to a physical state, b, is related 
through the orthogonal transformation, O: 

M b = J2°cbM a 

a 

= Y,°^e i5a e i4> -W a . (3.5) 

a 

In particular 

\M h \ 2 = E °lbWa + E 20 ab O pb cos(5 a + (t> a -8p- <f>p)W a W f) ,(3.G) 
\M T \ 2 = E W + E lOatOpb cos(<5 Q - (j> a - 5p + <f> fi )W a Wp.(3.t) 

a a>/3 

There is direct CP violation if 

\M b \ 2 - \Mj;\ 2 = -4 E O ab O pb sin(5 a - 5 P ) sin(0 a - <j> p )W a Wp (3.8) 

a>/3 
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does not vanish, which requires at least two strong channels, with different 
phase shifts. This could occur because two or more strong channels commu- 
nicate, or because there are two non-communicating channels that interfere 
(as in the case of K — > titt). In addition, there must be at least two distinct 
weak phases. Once both these conditions are met, in general there may be di- 
rect observable CP violation. If we sum over all final states, b, the difference 
in decay rates vanishes, as required by the CPT theorem. 

The strong-interaction matrix, S, has a number of invariant subspaces, 
that is, groups of states that don't get mixed with states from outside the 
group. Most obviously, isospin is a good quantum number so we can separate 
final states according to isospin. Similarly, parity and G-parity are good 
quantum numbers under strong interactions, and we can separate states by 
these, as well. Consider, for example, B° decays. The inr final states appear 
in both the 1 = and 1 = 2 subspaces. The 7T7T final states with J = have 
even parity and G-parity. The p-n final states are in different subspaces since 
they have odd parity and G-parity. The pp final states have even G-parity 
and 1 = and 1 = 2. However, they can have either even parity (s-wave 
and d-wave decays) or odd parity (p-wave decays). If we sum the partial 
decay rates over all final states with a given set of good quantum numbers, it 
must equal the rate of the decay of the conjugate particle into all final states 
with the CP-conjugate set of good quantum numbers. Here we must allow 
for mixing between states like DD with ttit, which although dynamically 
suppressed, is not forbidden through strong interactions. 

4 An Ansatz 

When the theory of final state interactions was first introduced 0, hadrons 
were considered elementary particles. Nonleptonic decays could then be 
viewed as a weak decay followed by hadronic final-state interactions. It 
was possible to imagine isolating the weak decay in the absence of strong 
interactions. This is no longer true for hadrons made of quarks, since the 
hadrons themselves are formed by the same forces responsible for the final- 
state interactions. 

To describe nonleptonic decays within the context of QCD we first inte- 
grate out the short-distance effects, generating an effective weak interaction. 
This introduces no ambiguity. However, it is still necessary to separate the 
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part of the long-distance interaction responsible for the formation of hadrons, 
from the remainder that causes the final-state interaction. 

Final-state interaction theory prescribes the form of the weak decay am- 
plitudes in the presence of strong interactions, but this, by itself, is not 
enough for us to draw useful conclusions. We adopt as our ansatz that the 
decay amplitude A4 a — e tSa e % ^ a W a is to be interpreted as a weak amplitude 
Ai a = e l ^ a W a with a strong final-state correction e tSa . The amplitude W a 
includes the physics of hadron formation, such as quark wave functions, as 
well as short-distance corrections to the weak interaction. Since the phase 
factor e l5a arises only from intermediate states that are on the mass shell, 
the off-mass-shell hadronic contributions are included implicitly in W a . With 
this interpretation, we can write the "bare" decay amplitude to the physical 
state b in the absence of final-state interactions as 

X b = J2°^ i(t>a W a , (4.1) 

a 

absorbing the weak phase in the definition of X a . Of course this is only 
heuristic: In the absence of strong interactions we wouldn't even have the 
hadronic final state b. Nonetheless, since we can write 

e^Wp = Y,°^ b , (4.2) 

b 

we can express the physical amplitude, M. in terms of the X amplitudes 

a a 

= Y,^baX a , (4.3) 

a 

where we have identified the "square root of the S-matrix" 

VS ba = J20 ab e^O aa . (4.4) 

a 

We shall take as our ansatz that we imagine the X b to be fixed, while we 
allow the strong interactions to vary, that is, while we allow the matrices O, 
and the phase shifts, S a , . . . to vary over some range. If the phase shifts are 
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all set to zero, there is no final state interaction and Aib = The opposite 
extreme allows the phase shifts to vary as widely as possible. 

The amplitude X b includes short-distance QCD corrections at energies 
much higher than the initial hadron mass, with the effective weak interac- 
tions improved by the renormalization group. In general, final-state inter- 
actions not only introduce phrases, but modify magnitudes as well. When 
we vary the phases randomly with our Ansatz, leaving the fixed, we are 
ignoring the necessary correlation between phases and magnitudes. This is 
an inevitable shortcoming of our statistical approach. 



5 Dyson's Random S-Matrices 

In 1962 Dyson proposed that in treating the distribution of resonance energy 
levels in a heavy nucleus, it would be better to consider random S-matrices 
rather than an ensemble of Hamiltonians. Because the S-matrices are uni- 
tary and thus bounded, they form a compact set and it is possible to define 
an invariant measure on them. The invariant measure provides the proper 
weighting to make the S-matrices random. Since the set of symmetric matri- 
ces do not form a group, finding the invariant measure required some careful 
analysis. Dyson was able to show that the proper weighting for a set of phase 
shifts (here we treat our "square root of the S-matrix" as the equivalent of 
Dyson's S-matrix) is 

V = [2 2JV 7r Ar r(l + \N)}- 1 [] \e iSi - e*H (5.1) 

i<j 

which has unit normalization: 

JHd5 k V = l. (5.2) 

k 

Altogether, the measure on the set of symmetric matrices is given by 

fi(d\^S) = n(dR) J] dS^V, (5.3) 

7 

where fi(dR) is the measure on the iV-dimensional orthogonal group. We 
shall never need the explicit form of fi(dR), but only the invariant nature of 
the measure. 
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If we postulate that the are fixed, while varying the strong interactions 
by integrating over the space of S-matrices (or rather their square roots), we 
can compute ensemble averages, at least in principle. 

Dyson's random S-matrices may be appropriate to decays in which final 
states acquire phases that may be large, that is, in the resonance region. The 
decays of D mesons might fit this description. When the final states are at 
higher energies, beyond the resonance region, we expect that the phases shifts 
may be small. This might describe B decays. Below we consider Dyson's 
completely random matrices and a number of alternatives, and calculate 
the observable quantities of interest for weak decays to see how final-state 
interactions manifest themselves. 



Dyson's construction of random S-matrices is a point of departure for phe- 
nomenological models. We can generalize by replacing V with some alterna- 
tive weight, W. The most obvious choice is the constant weight, Wconst = 



As noted in the introduction, the persistence of a correlation in the phase 
of different isospin amplitudes in B decays (the success of factorization) indi- 
cates that a completely randomized strong-interaction description cannot be 
successful. If the strong interaction effects become less important at higher 
energies, then a plausible description might be obtained by restricting the 
range of the eigenphase-shifts. We can model this behavior by selecting a 
function w and defining 



6 Models for Final State Interactions 



(2k) 



N 



(6.1) 
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with w(—S) = w(S) and 




'o 

2tt 




(6.2) 
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In the limit of small final-state interactions, 5 is confined to a symmetrical 
region around 5 = 0. At the other extreme, we can take w = 1/2tt and 
recover the results for W = W const - 

We can get a sense of these approximations by calculating ratio of the 
expectation values, averaged over the ensemble, of the elastic and total cross 
sections in the fixed zero angular momentum partial wave. Writing S = 
1 + 2iT, we have 

& elastic = ReT fe 2 b + ImT 6 2 b ^ ^ 

'total ImT 6 {, 

Since the S-matrix can be written as 

S ba = Y,0*be 2l5a O aa , (6.4) 

a 

we have 

ReS bb = ^O ab cos25 a O ab , 

a 

lmS bb = J20 ab sm25 a O ab , 

a 

(ReS bb ) 2 = J20 ab cos25 a O aa O pb cos25pOp b: 
(lmS bb ) 2 = J2°ab sin 28 a O aa O 0b sin 28 f3 O 0b . 

(6.5) 

Thus, for example, 

< {ReS bb f > = I n{dR) J] d6 y W £ O ab cos 28 a O aa O pb cos 25 p O pb . 

(6.6) 
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Integrals over the iV-dimensional orthogonal group of products of Os 
appear throughout our calculations. Results for up to eight Os are derived 
in Appendix [FJ. Using those results we find with W = W w 



< ReS bb > = < cos 25 >, 

< lmS bb > = <sin25>=0, 

< (ReS bb ) 2 > = — - — < cos 2 25 > + — — - < cos 25 > 2 , 
v ; N + 2 N + 2 

3 

< (lmS bb ) 2 > = j^—^ < sin 2 25 > . 



(6.7) 



Altogether, we find 



< ° elastlc > =< sin 2 5 > +— ^— < cos 2 5 > . (6.8) 

< Vtotal > N + 2 

For completely random phase shifts, < sin 2 5 >=< cos 2 5 >= 1/2, and the 
elasticity is equal to one-half, up to corrections of order 1/N. The observed 
ratio of elastic to total cross sections in np and Kp scattering at yfs ~ 6 GeV 
is about 15%. Of course, the scattering at these energies is due to many 
partial waves, but we should not be surprised to find a value of < sin 2 5 > 
as small as 0.1 to 0.2 at the B-decay energy. 

The average of the square of the decay matrix element over the ensemble 
of unitary (square-root of) S-matrices is given by 



< \M b \ 2 > 



f fi(dR) f n^wE (Vs ba x a ) £ (Vs bc x c 

J J a c 

= f p(dR) f l[d5 y Wj2T, c eb e l5a O aa ^20 pb e-^O^X a X* c . 

J J a,c a p 

(6.9) 

The simplest alternative is W = yV cons t, which has the virtue < cos 5 >= 
0. We proceed to calculate for this choice 

< \M b \ 2 > = f KdR) I II 7T 1 E E O ah e i5 «O aa ]T O&e-UeOpcXaX:. 

J 7 a,c a f) 

(6.10) 
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Only terms with a = (3 contribute, so 

< |7W 6 | 2 >= f nidR) E O ab O aa O ab O ac X a X* c . (6.11) 

a,a,c 

Using the results from Appendix [FJ, we find, for W = W cons t, 

< \M b \ 2 >= j^m 2 b \ + NW\), (6.12) 

where 

W\ = ^J2\X 2 \ (6.13) 

is the mean of the absolute squares of all the "bare" weak decay amplitudes. 

We now turn our attention to calculating the fluctuations away from this 
result, i.e., we consider 



<\M 2 b \ 2 > = J fx(dR) [HdSjW 



a,c,d,e a,/3,7,e 

xX a X* c X d X* e . (6.14) 
In evaluating this expression, we proceed by separating out categories with 



varying number of Greek indices are identical, as shown in Table ^TT| and 
discussed further in Appendix |Cj. 

If, as before, we set W = W cons t, we see that there are non-vanishing 
contributions only if a = 7 and 5 = e or if a = e and 7 = 5. There are two 
kinds of contributions: ones where all four Greek indices are the same, and 
ones where there are two distinct Greek indices. There are N of the former 
and 2N(N — 1) of the latter. Thus the sum over eigenchannels we need in 
calculating the average over the orthogonal group is 



n(dR) \N{O ab O aa O ab O ac O ab O ad O ab O ae ) 

+2N{N - 1) {O ab O aa O ah O ac Op h Op d Op h O Pe )}X a XlX d Xl 

(6.15) 
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Category 


Condition 


Terms 


Phase 


I 


a = 7 = 5 = e 


iV 


1 


II 


a = 7^(5 = e 


AT(AT-l) 


1 




a = e 7^ 7 = 5 


iV(iV- 1) 




III 


a = 5 7^ 7 = e 


N(N-1) 




IV 




N(N -1)(N -2)(N -3) 


g i(d"a — iy +<5<5 — tf e ) 


V 


a^7^(5 = e 


N(N-l)(N-2) 


gi(*a— A'-y ) 




5 7^ 7 7^ a = e 


N(N-l)(N-2) 


gi(5a-<5 7 ) 




5 = 7 7^ a 7^ e 


N(N-l)(N-2) 


gi(5 a — <5 £ ) 




a = 7 7^ 5 7^ e 


N(N — 1)(N — 2) 


gj(5 5 -<5 E ) 


VI 


a 7^ 5 7^ 7 = e 


iV(AT-l)(AT-2) 


gj(( ! > a +(3 7 — 25c ) 




a = 5 7^ 7 7^ e 


7V(7V- l)(JV-2) 


gi(25 Q -<5 7 -<5 (! ) 


VII 


a = 7 = 5 7^ e 


N(N-1) 


gi(#,5— <5 e ) 




a = 5 = e 7^ 7 


N(N-l) 


g«(5 a — <5 7 ) 




7 = 5 = e ^ a 


N(N-l) 


gi(5 a -<5 7 ) 




a = 7 = e 7^ 5 


N(N — 1) 


gj(55-<5 £ ) 



Table 6.1: Various categories of integrals over the iV-dimensional orthogonal 
group of products of eight Os. Since each Greek index can take on any of N 
values, the total number of terms is iV 4 . The integral over the phase depends 
on the choice of W. The result of the integral over the rotations is a tensor, 
which is contracted with X a X*X<iX*. 

Using the results of Appendix [FJ, and writing 

X~z = -J-£X 2 , 



N' 

X^ = ^EX* 2 . (6.16) 



for W = W, 



const i 



< \M 2 b \ 2 > 



(N + 1)(N + 2)(N + 4)(iV + 6) 

x 



2(iV + 3)(JV + 5)iV 2 (|X 2 |) 2 + 8N 2 {N + 5)|X 2 ||X fe 2 | 
-AN(N + 3)(X^X; 2 + X^X 2 ) + 8(N 2 + 7N + 9)|X 6 2 | 2 
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+N 2 (N + 9)X 2 X* 2 
To leading order in 1/N, the fluctuations are given by 



(6.17) 



< |A4 2 | 2 > - < |M 6 | 2 > 2 



2 

iV 1 



X 2 | + ^|X 2 | 



< \M b \ 2 > 2 



(6.18) 



It is not much harder to compute for the more general case W = W w . We 
find 



< |-M fe | 2 > 



N 



N + 2 
N 

N + 2 



-\X b \ 2 + |X 2 | + (|X b | 2 - |X 2 |) < cos<5 > 



< cos5 > 2 + — )\X b \ 2 + (1- < cos<5 > 2 )|X 2 | 



N 



(6.19) 



We recover the results for W = W CO nst by setting < cos 5 >= 0. In the limit 
of small phase shifts, we see that the effect of rescattering vanishes. 

A tedious calculation employing extensively the results of Appendix || 
yields the fluctuations in this result. To leading order in 1/N, we find, 



< \Ml\ 2 >-< \M b \ 2 > 



2^2 



IX 2 



1- < cos 5 > 2 ) 2 + (< cos 5 > 2 - < cos 25 >) 
+ |X fo | 2 |X 2 | 2 < cos5 > 2 —4 < cos5 > 4 +2 < cos 6> 2 < cos 25 > 



+\X, 



2 1 2 



N 



4 < cos 5 > 2 -8 < cos 5 > 4 +4 < cos 5 > 2 < cos 25 > 



(6.20) 



The fluctuations vanish, of course, in the limit of zero phase shift. In addition, 
the term varying as \X£\ is suppressed by 1/N. If the phase shifts are small, 



< \M b \ 2 >^ \X b \ 2 + <5 2 > |X 2 



(6.21) 



and the fluctuations are 
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< \M 2 b \ 2 > - < \M b \ 2 > 2 « 2<5 2 > 2 \X 2 \ +(<6 A >-<5 2 > 2 )(\X 2 \\X 2 \ + ^\X 2 \ 2 ). 

(6.22) 

We turn now to the calculation with the full Dyson weight, W = T> 



< \M b \ 2 >= / n(dR) / n^EE^ e40 ««E^" i,7 W* 

J J a,c a 7 

(6.23) 



We find after integrating over the Os: 



<i.M fe i 2 >= jnd6 k v 



N 



N + 2 



-\X b \ 2 + \X*\ + cos(5 1 - 5 2 )(\X 2 \ - |X2|) 



In @, Dyson shows how to calculate the n-point correlation function: 

N\ 

R n ($i, ...6 n 



. TT d5 k V. 
(N-n)\J £ 



(6.2 
(6.25) 



It is immediate that 

R 
Ri(S) 

In Appendix we show that 

/ ]JdS k Vcos(5 1 -5 2 ) 
J k 

Inserting this value, we find 



1. 

N 

2n' 



N 2 + 1 



(A^-1)(A^ 2 - 1) 



(6.26) 



(6.27) 



<\M b \ 2 >-- 



N 



N + 2 



N 



X b \ 2 + \X 2 \ 



N 2 + 1 



(N 2 - 1)(JV- 1) 



(I^|-I^ 2 l) 



(6.28) 

It is apparent that the inclusion of the Dyson factor induces only a small 
change, one which vanishes if the bare amplitude has the same magnitude as 
the average amplitude, and which is suppressed by one power of 1/N. 
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7 Coherence of Interfering Amplitudes 



Coherence of interfering amplitudes is a measure of final state interactions. 
We can address this question quite simply using the weight W w . Suppose 
we have two interfering amplitudes from distinct channels (e.g., I = 1/2 and 
1 = 3/2) so that the full amplitude is 

M b + Mv=J2 O ab e iS "O aa X a + 0' a , b ,e^O' a , a ,X a , (7.1) 

a,oi a', a' 

When we average over the ensemble of S-matrices we average over the two 
sectors containing b and b', which are completely independent. We thus find 



< \M b + M b >\ 2 > 



N 



N + 2 
N' 



< cos 5 > 2 +-=)\Xg\ + (1- < cos 5 > 2 )\X 2 



+ 



N> + 2 



< cos 5' > 2 +—)\X 2 \ + (1- < cos 5' > 2 )\X' 2 \ 



+2ReX b Xt, < cos 5 >< cos<5' > . 



(7.2) 



If the phase shifts are small, coherence is maintained as is seen from 
rewriting this as 



< \M b + M b >\ 2 > 

= | < COS (5 > X b + < cos 5' > X v \ 2 
iV(l-<cos^> 2 ) / 2 — 
+ N~T2 U |Xfel + |X 1 

iV ' a - <cos5 ' >2) f hx 2 \ +W \ 



N' + 2 VA^ 



(7.3) 



When the phase shifts vary over a broad range so that < cos 5 >~ 0, the 
right hand side of Eq. ( |7.3| ) gives simply the incoherent sum: 

< \M b + M b >\ 2 >^< \M b \ 2 > + < \M b >\ 2 >; <cos5>^0, (7.4) 
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and there is no observed interference. In contrast, when the final state inter- 
actions are small and < 5 2 >^ 0, 



< \M b + My\ 2 > — > I < cos<5 > X b + < cos<5' > X b ,\ 2 

+ <5 2 > \X T \ + < 5' 2 > \X^\ 

+0(1/N); < 5 2 >^ 0. (7.5) 

Numerically, the incoherence introduced by final state interactions in this 
case turns out to be surprisingly small. It is the feature that makes the 
factorization hypothesis work in the B decays, with their small final state 
interactions, whereas factorization fails in the D decays, where the final state 
interactions are substantial. A more quantitative treatment of these issues 



is given in Section 10 



8 Direct CP Violation 

As discussed in Section 2, direct CP violation is observable when the X a are 
complex, provided there are final state interactions and more than one weak 
phase. The rates for a decay and its CP conjugate are given by 

<|-M 6 | 2 > = f fi(dR) f]Jd5 k WY / T, -be iSa O aa Y / O lb e- l5 -'0 lc X a X:, 

J J a,c a 7 

a,c a 7 

(8.1) 

These have equal expectations provided there is no a priori distinction 
between the channels, i.e. so long as the weight W(5i, S 2 , ■ ■ ■ <5jv) is symmetric 
under the interchange of any two Ss. This is true of all the models we consider 
and is implicit in our statistical approach in general. Of course, any particular 
choice of S-matrix will introduce direct CP violation, but the effect averaged 
over any ensemble symmetric in the 5s will vanish. In particular, we can 
write 

a,c a 7 

(8.2) 
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Consider next 



< (\MI\ - \Ml\f > 

= J nidR) J ]Jdd k W £ o ab o aa o, b o, c o 5b o Sd o tb o ee 

a,b,c,d o,/3,7,e 

x 2 e <(*«-*r)( e <fo-*«) _ e - i{Ss - Sc) )X a X*X d X*. (8.3) 

The effect of the second term in the bracket is to change the phase in a 
particular way. Referring to Table |TT] and Appendix |Cj, we see that the 
phase change interchanges the second instance in Category II with Category 
III, and the second and third instances of Category V to Category VI. To 
leading order in N, then, the result is, for W = W w , 



< (\Ml\ - \M 2 b \) > = 2(|X 2 |-X 2 X* 2 )(l-<cos25> 2 ) 

= 4* E |^a| 2 |X c | 2 sin 2 (0 a - C )(1- < cos 25 > 2 ). 

a,c 

(8.4) 

This is to be compared to Eq.( |6.19| ) taken in the same limit, 

< \M b \ 2 >=]X 2 ~\+ <cosd> 2 (|X 2 | -X 2 ). (8.5) 

With large final state interactions so < cos 5 >= 0, the relative magnitude 
of CP violation is 



< (\M 2 b \ - \M 2 \f > E a , c |X a | 2 |X c | 2 sur 



2 ^ , T ^2 



< (J-Mfl + \Ml\) > iV 2 X 2 

(8.6) 

When the phase shifts are small, < 5 2 ><^C 1, the relative magnitude is 



<( 




- \Ml\ 


) 2 > 4<6 2 >Ea,c\Xa\ 


2 \Xt 


: | 2 sin 2 (0 a -0 c ) 


<( 




+ \Ml\ 


f> " mx 2 b \- 


2 
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which, naturally, is suppressed by the small strong phases, as anticipated in 
Eq. §3). 



9 Relation to Quark-level Analysis 

The short- distance portion of final state interactions has been computed in 
perturbative QCD. The result has been applied to CP violation phenomena 
in B decays using some models for the soft strong interactions ||, |9|, |lOfl . 
Our statistical approach is complementary to this QCD based model. We 
formulate final state interactions solely in terms of states of hadrons, filling in 
for our ignorance of detailed interactions with the assumption of statistical 
randomness. Rather than trying to predict the precise value of a process 
we attempt to estimate an ensemble average, the ensemble being a range of 
strong interactions, as characterized by their S-matrices. 

In the quark-level approach, the effective weak interaction Hamiltonian 
is written 

n = J2^)o l , (9.8) 

i 

where Ci{fj) incorporates the short-distance QCD loop corrections, and the 
Oi are four-quark operators that carry weak CP-violation phases. The co- 
efficients Cj(/i) are computed by the renormalization group method with fj, 2 
kept in the spacelike region. Therefore the Ci(fi) are real and thus provide 
no final-state-interaction phases. If the matrix elements of the operators Oi 
are computed allowing one-loop corrections, imaginary parts can arise. The 
one-loop corrections of the matrix elements combine with corrections to the 
coefficient functions so that the result may be written 

H = J2cf f (riOt ee , (9.9) 

i 

where now the q do carry strong interaction phases. Though the one-loop 
corrections to the matrix elements produce strong interaction phases for all 
interactions Oi, Ref. [10] argues that the dominant phases arise from matrix 



elements of penguin-type operators, the source identified in Ref. . 

This procedure generates final-state interactions from short- distance phe- 
nomena. The formation of the hadrons itself introduces no phases, since it 
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is carried out using models [|Tl|] based on factorization. The accuracy of this 
perturbative calculation of strong phases is yet to be tested, but it has been 
advocated on the basis of asymptotic freedom 0. 

In the quark-based approach, the strong phase of < b\H e ff\B > arises 
from the effective coefficient functions of the various four-quark operators. 
Because different final states receive different contributions from tree opera- 
tors, their phases, which are generated by the coefficients of these operators, 
will vary. 

A rough comparison can be made with our approach. We might imagine 
the effective Hamilton of Eq. ( |9.8| ) to give our amplitude without final-state 
interactions, i.e. Eq. Q4.1|) . In place of the perturbative calculation of hadronic 
rescattering, our procedure introduces the eigenphases of the S-matrix, which 
are, of course, unknown. This results in a loss of explicit predictive power, 
but provides a convenient formalism for discussing final-state interactions. 



10 Comparison with Experiment 

We have argued that a purely random S-matrix, with the phase shift taking 
on all possible values with equal likelihood would work best at energies where 
the strong interactions produced resonances, while above those energies the 
phase shifts would be confined to relatively small values. Thus in the decay of 
D mesons, with mass a bit below 2 GeV, the random model might apply, while 
for B meson decays, the final state interactions should be less substantial. 
We consider D and B decays in light of these expectations. 



10.1 D decays 

The development of the preceding Sections images all decays to be into two- 
body final states. Even broadening this to include quasi-two-body states 
cannot account for all decays. How then are we to determine N, the number 
of discrete final states? We work empirically, inferring an effective channel 
number that acts like N. Fortunately, our predictions do not depend crucially 
on the value of N, only on its being large. 

We estimate N by looking at the observed branching ratios. The fully 
random model predicts that the squares of the matrix elements should cluster 
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around the mean 

<|7W f) | 2 >^p^ (10.10) 

with the standard deviation equal to the mean itself. Continuing to neglect 
phase space corrections, we can relate the mean branching ratio to the inclu- 
sive branching ratio to the category considered. Thus the effective number 
of decay channels for non-leptonic D decays is roughly 

N _ BR ( inClUSiVe )n 0n -Ze P tom C ^ Q ^ 

B Rnon—leptonic 

Here we ignore the Cabibbo-suppressed decays. We must, however, separate 
the + and 0~ states since they are not mixed by the strong interactions. 
The relevant branching ratios taken from the Review of Particle Properties 
[§] are listed in Tables ffO and [TO. 



D+ 


BR in % 


D+ 


BR in % 


J p = + final state 




J p = 0~ final state 




17° + 


2.74 ±0.29 


1P U + 
K 


6.6 ±2.5 


IF + 


8.1 ±1.7 


K aj 


< 0.3 


K p + 


2.1 ± 1.4 




1.92 ±0.19 


K 1 (1270)°-k + 


< 0.7 


Z(1410)*°7T + 


< 0.7 


Fi(1400)°7T + 


5.0 ± 1.3 


Z (1430)* 7T + 


3.7 ±0.4 






Z(1680)*°7T + 


1.45 ±0.31 



Table 10.1: Branching ratios of D + into states with J p = + and J p = 0~ 
taken from Ref. ||. Final state that can be either + or 0~ are also included 
under + . The mean branching ratio for + is 3.6 ± 0.9%. while for 0~ it is 
2.3 ±0.6%. 



The inclusive non-leptonic branching ratios for neutral and charged D 
mesons are about 0.8 and 0.6 respectively. If we ascribe half of each to- 
tal to - and half to + , then for the charged D's, we find N 3 / 2 ~ 8 for 
the parity-even states and N 3 / 2 ~ 13 for the parity-odd states. In both 
cases, the final states have / = 3/2. The neutral D decays have the added 
complexity that they involve both 1 = 1/2 and / = 3/2. Thus we infer 
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D° 

J p = + final state 


BR in % 


D° 

J p = 0~ final state 


BR in % 


K-7T + 


3.83 ±0.12 


K-p + 


10.8 ± 1.0 




2.11 ±0.21 




1.20 ±0.17 




0.7 ±0.1 




2.1 ±0.4 


i\ // 


1 70 + 26 


if f'(QRQ) 


57 + 1 6 




70J.1 1 
( .O ± 1.1 




U.oD ± U.1U 


17° 


< 1.9 


K / 2 (1270) 


r\ A "i 1 /a -1 r" 

0.41 ± 0.15 


Try 


1.47 ±0.33 


K-a 2 (1320) + 


< 0.2 




6.0 ±2.4 




5.0 ±0.4 


F x (1270)-7r + 


1.06 ±0.29 




3.1 ±0.4 


iri(1400)-7r+ 


< 1.2 


K(1410)*-7r + 


< 1.2 


ir 1 (1400)V ) 


< 3.7 


K (1430)*-7r + 


1.24 ±0.26 




1.1 ±0.4 


ir 2 (1430)*-7r+ 


< 0.8 






K(1410)*°7r° 


< 0.4 






K*°7] 


1.9 ±0.5 






K*°rf 


< 0.11 



Table 10.2: Branching ratios of D° into states with J p = + and J p = 0~ 
taken from Ref. ||. Final states that can be either + or 0~ are also included 
under J p = + . The mean branching ratio for J p = + is 2.1 ± 0.4%, while 
for J p = Q- it is 1.7 ± 0.2. 



N3/2 + N 1/2 « 0.4/0.02 « 20. Roughly, then, N 1/2 and N 3/2 are each about 
10, for both the parity-even and parity-odd sectors. In Fig. [I] we show 
the distribution of branching ratios for D + and D°. In Fig. || the data 
are combined and compared with an exponential distribution of the form 
exp(— |-Mfc|/|< >)• This form has the virtue of satisfying Eq. ( |6.18|) , 
relating the standard deviation of the distribution of |-M;,| 2 to its mean. The 
agreement with this empirical form is generally satisfactory. If we compute 
the standard deviation of the branching ratios, combining data from all de- 
cay channels of neutral and charged Ds and ignore phase space corrections, 
we find 
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A\Ml\ 



;i.2±0.2) < \Ml\ > 



in agreement with our expectation A|.M 2 | 



=< \Ml\ 



>. 
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(b) 



n 



6 8 10 12 



Branching Ratio 



Branching Ratio 



Figure 1: (a) Decay histograms for D + non-leptonic decays, combining both 
J p = + and 0~. The number of modes have a branching ratio between 0% 
and 1%, between 1% and 2%, etc. is displayed. (b)Decay histograms for D° 
non-leptonic decays, combining both J p = + and 0~. 



10.2 B decays 

Because of the larger phase space, the number of channels available for B 
decays should be substantially larger than for D decays. This view is sup- 
ported by the absence of any dominant exclusive B-decay modes. We expect 
the large N limit to be very good in B decays, but because of the greater 
energy available for decays, weaker final state interactions are anticipated. 
We address this by restricting the range of S. 

The value of sin 2 5 is equal to the ratio of the elastic to the total cross 
section in a fixed partial wave by Eq. fl6.8|) , which to order (1/N) gives 
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2 4 6 8 10 12 
Branching Ratio (%) 

Figure 2: Decay histogram for D + and D° non-leptonic decays, combining 
both J p = + and 0~. The curve has an exponential fall-off with the standard 
deviation equal to the mean, as predicted by the purely random phase shift 
model. 



sin 2 5 = ^± ( 10 .13) 

<7 tot 

Of course experiments are performed with beams, and therefore measure the 
sum of many partial waves. If we make the crude approximation of assuming 
the same elasticity for all partial waves, we may estimate < sin 2 5 0.2, 
on the basis of typical hadronic cross sections evaluated for a cm. energy 
around 5 GeV. 

Consider two different distributions of 5, one flat between — 5 max and 

& max-, 

dP 1 



d5 25 max ' 
and the other one Gaussian 



-8 max < 6 < 5 max (10.14) 



dP 1 

^ = ^exp(-5 2 /5 2 ). (10-15) 
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If we fix 5 m ax = 1-66 and 5o = 0.715, then both give < sin 2 8 >= 0.2. The flat 
distribution gives < cos 5 >= 0.89, while the Gaussian gives < cos 5 >= 0.88. 

We can apply such models to estimate the breakdown of coherence (which 
is essential to factorization) caused by final state interactions. We can write 
Eq. O as 



< \M b + My\ 2 > 

=< cos 5 > 2 \X b + X b ,\ 2 + (1- < cos 5 > 2 ) p^ + ^l] ,(10.16) 

where we have dropped terms of order 1/N and taken the distributions of 
5 and 5' to be identical. Factorization leads to the prediction that some 
branching ratios are much smaller than their isospin relatives, e.g., BR(B° — > 
D°p°) <C BR(B° — > D~p + ). Thus for one isospin variant we have X b + 
Xy m 0, while for another variant we expect X b + Xy to be typical of modes 
generally, i.e. roughly \J\X 2 \ + \J\X' 2 \. 

If we suppose complete cancellation for one mode without final state 
interactions, then the ratio of the decay rate into that mode divided by the 
rate into the isospin variant can be written 

1 " <C0S5>2 (10.17) 



1- < cos<5 > 2 + < cos5 > 2 j=Sffi' 

where the amplitudes X b and Xy are pure isospin amplitudes for the unsup- 
pressed mode. If all amplitudes have about the same size, this ratio is simply 
(using the approximate value of < cos 5 > above) 

1— < cos 5 > 2 

-^0.13, 10.18 

1+<cos5> 2 

whereas if the modes in question have branching ratios twice as as large as 
typical modes, the ratio is 

1 " <C0S5>2 0.07. (10.19) 



1 + 3 < cos 5 > 2 



It is not possible in the case of B decays to estimate the relation between 
any particular amplitude and the mean of the amplitudes as in the case of 
D decays. The complete mixing of the amplitudes by final state interactions 
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in the case of D decays makes each one a better representative of the whole 
group. 

The current experimental values for the B — > Dp system are || 

BR(B° -> D°p°) < 5.5 x 10" 4 (10.20) 

and 

BR(B° -> D-p + ) = (7.8 ± 1.3) x 10~ 3 . (10.21) 

Final state interactions are not the only contribution to decays like B° — > 
D°p°. The small value of the ratio BR(B° -> Tfp°)/BR(B° -> £>-p+) in- 



dicated by Eqs. ( pLQ.20| ) - ( |10.21| ) suggests that the final state interactions, 
as characterized by < sin 2 5 > are probably not too large. Of course, even 
if there are large final-state interactions, by chance the phases for ampli- 
tudes contributing to a single final state might end up coherent just as if 
there were no final-state interactions. This may be a case for the processes 
D° — > K~p + and D° — ► K°p°, whose large ratio agrees with the factorization 
prediction. However, the failure of factorization in the D° — > K*n decays 
is enough to demonstrate that final-state interactions are important in D 
decays, supposing that there is some validity to the factorization approach. 
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A Final-State Interaction Theorem 



We demonstrate here the final-state interaction theorem in the form appro- 
priate to our considerations of weak decays. We begin with an iV-channel 
strong interaction S-matrix, S, and then add a single channel that is coupled 
to the others only by weak interactions. We write the expanded S-matrix as 

s '={u C ' T s)- 

where U is an iV-element column vector, and U T is its transpose. We have 
assumed T invariance so that the full S-matrix is symmetric. Unitarity tells 
us that 



1 + U T U* = 1, (A.2) 
U + SU* = 0, (A.3) 
SS ] + UU ] = I. (A.4) 

Since U is a weak amplitude and thus small, we drop the U 2 terms and 
are left with a single relation, U = —SU*. Because S is not just unitary, but 
symmetric, it can be diagonalized by an orthogonal (rotation) matrix: 

S diag = OSO T , (A.5) 

where Sdiag has diagonal elements e 2tSa ;a = l,...N. Now our equation is 
simply 



OU = -S dtag OU*. (A.6) 



If we set V = OU, then 



and 



V a = -e 2lS «V:, (A.7) 



V a = ie i5a R a , (A.8) 

where R is an arbitrary real number. The subscript a labels the individual 
eigenchannels of the S matrix. In particular 



28 



U b = iJ20a b e iSa R, 



(A.9) 



a 



where we have used Latin subscripts to indicate physical states, as opposed 
to the eigenchannels labeled in Greek. 

Thus, in the case of no CP violation (so time reversal invariance is true) , 
final state interactions provide phases as follows. Each eigenchannel of the 
S-matrix gets some phase. The physical states are combinations of the eigen- 
channels and pick up corresponding phases. 

How is this changed by CP violation? We no longer can assume the full 
S-matrix is symmetric, but we can assume the strong interaction S-matrix is 
symmetric. We write in this instance 



The idea is that the weak amplitude has weak phases, which come in with 
one sign in LL and the other in LL. With this convention, LL gives the matrix 
elements for the decay of the single-particle state, while LL gives the matrix 
elements for the multiparticle states to reform the single-particle state, or for 
the time-reversal of the single-particle decay. Since CP is equivalent to T, it 
is clear that LL gives the amplitudes for the CP-conjugate decays. Now we 
find from unitarity 




(A.IO) 



LL + SU*_ = 0, 
LL + SU1 = 0, 



(A 
(A 
(A 



11) 
12) 

13) 



or equivalently 



U+ + U. 
LL - U. 



-s(u + + u.y, 

S(U+ - LL)*. 



(A 
(A 
(A 



14) 
15) 
16) 



Proceeding analogously to our earlier development, 
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OU ± , (A.17) 
-e 2i5a (V + + V-)* a , (A. 18) 

e 2iSa (V + -V_)* a , (A. 19) 

(A.20) 

and therefore 

(V+ + V-) a = 2ie iS "R a , (A.21) 
(V+-V-) a = -2e i5 «I a , (A.22) 

(A.23) 

where R a , I a are real numbers. Altogether 



V± = 
{V+ + V.) a = 
(V + -V) a = 



V ±a = ie iSa (R a ±U a ) =ie iSa e ±i<l,a W a , (A.24) 
U ±b = iY J O ab e i5 -{R a ±U a )=iY J O ab e i5 -e ±i ^W a . (A.25) 

a a 

When an approximate selection rule or dynamical suppression mechanism 
strongly hinders transitions between two groups of states, it is sensible to 
treat transitions between the two groups as a perturbation. For example, 
transitions among the states if}7nr, DD, and tttt are all allowed by strong 
interactions, but inhibited dynamically. We can represent such a state of 
affairs by an S-matrix that is nearly block diagonal. Consider a circumstance 
in which two different groups having Ni and N 2 channels are connected only 
slightly by the strong interactions, so the full S-matrix is 



S' 



( i XJl \ 

Z7i SV> e , (A.26) 
V U 2 e T J 

where we have assumed CP invariance for the present. We can diagonalize 
SW and ^ : 

= OWSVOW; S®, = O^SWK (A.27) 
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If we ignore terms of order e 2 , the diagonal elements of S^ iag and S^ iag are 

of the form e lS "\ a = l,...Ni and e s P , (3 = 1, ...N2. In the diagonalizing 
basis, as a result of unitarity the mixing matrices satisfy 



7 

X! e 7/3 e 7/3' = 0; P^P', 



7 



and they can be written in the form 



Cap = l^ape 



where e a/ 3 is real. 

Writing Vi = 0®Ui, we find 



2Re(e-^V«)+^e^V / S 2) * = 0, 
2Re(e-<V /3 (2) )+^e^ 1 Vi 1 )* = 0. 



The solution is 







(A.28) 



«( s » +i D, (A.29) 



(A.30) 



(A.31) 



where W^ 1 ** and W^ 2 ) are real. Unitarity provides no information about the 
real part of the amplitude (which here is proceeded by the i that connects the 
S matrix with the invariant amplitude). It is proper, however, to consider 
to have contributions of order e that arise from the mixing perturbation. 
The reaction rate will thus receive corrections of order e and of order e 2 from 
the real, undetermined part, and corrections of order e 2 from the imaginary 
part that are explicitly determined. 
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If CP violation is included, the amplitudes W^' acquire weak phases that 
change sign under CP. The result is 



P 



(A.32) 



This can be made more concrete by considering a salient example. The 
decay B — > 7r + 7r~ is important for the measurement of sin 2a, where a is 
one of the angles of the unitarity triangle. The decay receives a tree-level 
contribution from the quark-level process b — > uud, which has an amplitude 
proportional to V u bV* d . Were that the only contribution, there could be no 
direct CP violation, which requires at least two weak phases. On the other 
hand, CP violation could be observed through B°B mixing, and that CP 
violation would be especially easy to interpret because it would be unaffected 
by final state interactions. However, the irn final state is connected by strong 
interactions to DD. The decay B — > DD proceeds at tree-level through 
b — > ccd, with an amplitude proportional to V c bV* d . Since V u b/V c b ~ 0.08, 
while V c d/V U( i ~ 0.2, the overall process B — > DD — > tttt may be important. 
This is precisely the process known colloquially as the "penguin." Of course, 
the state DD represents the totality of other similar states, DDn, DD7nr, 
etc. 

In the perturbative quark-level calculation, the penguin diagram con- 
tributes a part that is real (aside from the Kobayashi-Maskawa factor), anal- 
ogous to the implicit contribution to W^ 1 ' above, and an imaginary part[[3| 
analogous to the explicit mixing displayed in Eq. (|A.32p . In the quark- 



level calculation, this absorptive contribution from the penguin diagram is 
the only source of final state interactions. In our formalism, the final state 
interactions are manifested the various phases of the eigenchannels. 



B Integrals over the Orthogonal Group 

We demonstrate here the computation of various invariant tensors obtained 
by integrating over the N dimensional orthogonal group. 
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If we consider the tensor 

J n(dR)O aa O pb , (B.l) 

we see that 

This follows from the invariance of the measure, fi(dR): we can shift the 
integration variable from O to O'O. A similar manipulation shows that 

Taai3bO' ac O' bd = T ac/3d . (B.3) 

In other words, the tensor T is invariant. We assume below that N > 8. 
This eliminates, for integrals of up to eight Os, the possibility of invariants 
formed from the determinant and means that every tensor must be made 
from products of Kronecker 5s that combine either first indices (which are 
indicated by Greek letters) or second indices (which are indicated by Roman 
letters). 

For example, 

J ii{dR)O aa Op h (B.4) 

must be proportional to 5 a pSab- Contracting the integrand with 5 a b gives 
unity, so 

J n{dR)O aa O pb = ^5 aP 5 ab . (B.5) 

Consider next 

J ix{dR)O aa Op h lc 5d , (B.6) 

which must be constructed from terms like 5 a p5^s5 a b^cd^ S ai 5psS a b5 C( i, etc. 
There are precisely nine such terms since there are 4!/(2!2!2!) = 3 ways to 
pick the Greek-index tensor and an equal number of ways to pick the Roman- 
index tensor. With so few terms it is easy to perform enough contractions 
to show that 

J n(dR)O aa Op b lc O Sd 
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— (SapS-ysSaeSbd + &*/3&y5 ^ad&c + S ay 5psSabS c d 

In particular, for ct 7^ 7. 

— jy^jy _ _|_ 2) ^^ a l a 2^ a 3 a 4 — ^aia 3 ^a 2 a4 — ^aia4^a 2 a 3 ] (-B.8) 

A straightforward extension of this technique to calculate the integral of eight 
Os leads to consideration of [8!/(2!2!2!2!4!)] 2 = 105 2 = 11,025 terms, so a 
different strategy is required. 

Consider the class of tensors formed when all the Greek indices are iden- 
tical. Now consider 

f 4 

/ n(dR) O aai = A^{5 aia2 5 a . iaA + S aias 5 a2a4 

+ S aia J a3a2 ). (B.9) 

J 1=1 
Multiply by 5 a3 a 4 and find 

/ fi(dR) n o aai = h aia2 
j i=i jv 

= A 4 (N + 2)5 aia2 , (B.10) 

so that A A = l/[JV(JV + 2)]. Using similar definitions and procedures we find 
analogously 

6 



U(dR)l[O aai = A 6 Yl 

i=l pairings 

1 

6 ~ N(N + 2)(N + 4)' 

8 

j n(dR)Y[O aai = A 8 Sa ta:i 5 akai 5 aman 5 apaq , 

i=l pairings 

1 



iV(iV + 2)(iV + 4)(iV + 6)' 
34 



(B.H) 



where the product over pairs has 105 = 8!/(2!2!2!2!4!) terms. 

Consider the integral of the product of six Os. The case of six identical 
Greek indices is computed above. There are two other possibilities: four 
identical indices and one distinct identical pair or three distinct pairs. 

Let us write (ij) for S aia .. Then we can represent the integral of six Os, 
with four identical Greek indices as 

J fl(dR)O aai O aa2 O aa . a O aa4 0/3a 5 Of3 a(i 

= A 42 {(12)(34)(56) + (13) (24) (56) + (14)(23)(56)} 

+B 42 {(12)(35)(46) + (12)(36)(45) + (13)(25)(46) + (13)(26)(45) 

+(14)(25)(36) + (14)(26)(35) + (15)(23)(46) + (15)(24)(36) 
+(15)(34)(26) + (16)(23)(45) + (16)(24)(35) + (16)(34)(25)} , 

(B.12) 

where the twelve B 42 terms mix 5 and 6 with 1,2,3,4. Contracting the 
integral over six O's with (56) gives, 



A 4 [(12)(34) + (13)(24) + (14)(23)] 

= At 2 iV[(12)(34) + (13)(24) + (14)(23)] 

+4 J B 42 [(12)(34) + (13) (24) + (14) (23)], (B.13) 

so 

NWT2)= NA " + iB "- (B14) 
Contracting next over (34) gives from Eq.( B.7| ), 

' -[(iV+l)(12)(56)-(15)(26)-(16)(24)] 



N(N — 1)(N + 2) 

= A 42 (N + 2)(12)(56) + 5 42 {2(12)(56) + (4 + iV)[(15)(26) + (16)(25)]}. 

(B.15) 



Thus 

N+l 



(N + 2)A 42 + 2B 42 , 



N(N - l)(N + 2) 

' (N + A)B 42 , (B.16) 



N(N -1)(N + 2) 
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so we have 

iV + 3 1 



.4 



12 



N-l N(N + 2)(N + 4)' 



Bi2 (N- 1) ' N(N + 2)(N + 4)' ^ B ' 17 ^ 
We proceed similarly when there are three distinct pairs of Greek indices: 

J fi(dR)OaaiO a a 2 @ 'f3a 3 '/3o 4 Ojas O yafi 

= A 222 {(12)(34)(56)} 

+B 222 {(12)(35)(46) + (12)(36)(45) + (15)(34)(26) 

+ (16)(34)(25) + (13)(24)(56) + (14)(23)(56)} 

+C 222 {(13)(25)(46) + (13)(26)(45) + (14)(25)(36) 
+ (14)(26)(35) + (15)(23)(46) + (15)(24)(36) 
+ (16)(23)(45) + (16)(24)(35)}. 

Contracting with (56) and comparing with Eq.( |B.7| ) we find from the 
coefficients of (12) (34) and (13) (24) 

N + l 

NA222 + 4B222 



N(N -1)(N + 2)' 

NB222 + 4C 222 = 



iV(JV-l)(JV + 2)' 

Contracting with (45) gives 

A 222 + {N + 1)B 2 22 + 2C 22 2 = 0, 
2B222 + (N + 2)C 222 = 0, 

so that 

N 2 + 3N -2 



(B.18) 



(B.19) 



.4 



222 



N(N - 1)(JV - 2)(N + 2)(N + 4) 
N + 2 

-D222 — 



N(N- 1)(N - 2)(N + 2)(N + A)' 
36 



° 222 N(N - 1)(N-2)(N + 2)(N + A)' 



(B.20) 



For the integral of eight Os, we have five cases, indicated symbolically by 
8, 62, 44, 422, and 2222, where the numbers represent the number of identical 
Greek indices. The '8' case is done above. For '62' we write 



J Li {dR^ (D aai O aa2 O aa3 O aai O aa5 Oaatj @/3a7 0/3as 

= A 62 {(78) [(12)(34)(56) + (12)(35)(46) + (12)(36)(45) 
+ (13)(24)(56) + (13)(25)(46) + (13)(26)(45) 
+ (14)(23)(56) + (14)(25)(36) + (14)(26)(35) 
+ (15)(23)(46) + (15)(24)(36) + (15)(26)(34) 
+ (16)(23)(45) + (16)(24)(35) + (16)(25)(34)]} 

+B 62 {(17)(28) [(34) (56) + (35)(46) + (36)(45)] 

+ (17)(38) [(24)(56) + (25)(46) + (26)(45)] ...} . 

(B.21) 

First contract with (78). Then identify the coefficient of (12) (34) (56): 

+ = N(N + l ){N + i y (B.22) 
Next, contract with (56) and identify the coefficient of (12) (34) (78). 



(N + 4)A 62 + 2B 62 = A 42 = —^-- — ^77 C • (B.23) 

V ; 62 62 42 N-l N(N + 2)(N + 4) V ; 

It follows that 

N + 5 



A 



62 



N(N - 1)(N + 2)(N + 4)(N + 6)' 
1 

r> 

62 N(N -l)(N + 2){N + 4){N + 6) 



[B.24) 
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Now let us deal with the 2222 computation. 



J fl(dR) O aai O aa2 Op a „ 0/3a 4 Oya 5 C 7 a 6 0$a 7 Osag 

= A 2222 {(12)(34)(56)(78)} 

+5 2222 {(12) [(35) (47) (68) + (35)(48)(67) + (36)(57)(48) + (36)(58)(47) 
+ (37) (46) (58) + (37) (45) (68) + (38) (46) (57) + (38) (45) (67)] 
+(34) [(15) (27) (68) + (15)(28)(67) + (16)(57)(28) + (16)(58)(27) 

+ (17)(26)(58) + (17)(25)(68) + (18)(26)(57) + (18)(25)(67)] . . .} 
+C2222 {(12)(34)[(57)(68) + (58)(67)] + (12) (56) [(37) (48) + (38)(47)] 
+(12) (78) [(53) (64) + (54)(63)] + (34) (56) [(17) (28) + (18)(27)] 
+ (34)(78)[(15)(26) + (16)(25)] + (56) (78) [(13) (24) + (14)(23)]} 
+^2222 {[(13)(24) + (14) (23)] [(57) (68) + (58)(67)] 
+ [(15)(26) + (16)(25)][(37)(48) + (38)(47)] 
+ [(17)(28) + (18)(27)][(53)(64) + (54)(63)]} 
+£2222 {(13) (25) (47) (68) + (13) (25) (48) (67) + (13) (26) (47) (58) + (13) (26) (48) (57) 
+(13) (27) (45) (68) + (13) (27) (46) (58) + (13) (28) (45) (67) + (13) (28) (46) (57) 
+(14) (25) (37) (68) + (14) (25) (38) (67) + (14) (26) (37) (58) + (14) (26) (38) (57) 
+ (14) (27) (35) (68) + (14) (27) (36) (58) + (14) (28) (35) (67) + (14) (28) (36) (57) . . .} . 

(B.25) 

First contract with (78). From the coefficients of (12) (34) (56), (12) (35) (46), 
and (13) (25) (46), respectively, we find 



NA2222 


+ 6C 2222 


= A222 


4-B 2222 + NC2222 


+ 2Z) 2222 


— B222 


NB2222 


+ QE2222 


= C222 



(B.26) 

Contracting, instead, with (67), we find 

^2222 + 45 2222 + (N + 1)C 2222 = 0, 
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(TV + 2)Z?2222 + 2C2222 + 2i?2222 = 0, 
273 2 222 + £>2222 + (TV + 3)i?2222 = 0. 

(B.27) 



The solution is 

(TV - 2) (TV 3 + 97V 2 + 197V + 3) 



^2222 = 

-B2222 = 

C2222 = 

-D2222 = 

-^2222 = 



N(N - 1)(7V + 1)(7V - 2)(/V + 2)(7V - 3)(/V + 4)(7V + 6) ' 

27V(7V + 4) 

7V(7V - 1)(7V + 1)(7V - 2)(7V + 2)(7V - 3)(7V + 4)(7V + 6) ' 

TV 3 + 67V 2 + 37V - 6 
~7V(7V - 1)(7V + 1)(7V - 2)(7V + 2)(7V - 3)(iV + 4)(7V + 6) ' 

TV 2 + 57V + 18 

7V(7V -1)(N + 1)(7V -2)(N + 2)(7V - 3)(7V + 4)(7V + 6) ' 

57V + 6 

~7V(7V - 1)(7V + 1)(7V - 2)(7V + 2)(7V - 3)(7V + 4)(7V + 6) ' 

(B.28) 



Turning next to 422, we write 

J fi(dR) O aai O aa2 O aa . A O aa4 OfS a5 Op a() O ia7 O iag 
= A 422 {[(12)(34) + (13)(24) + (14)(23)] (56)(78)} 

+5422 {[(12) (34) + (13)(24) + (14)(23)] [(57)(68) + (58)(67)]} 
+C422 {(78) [(15)(26)(34) + (16)(25)(34) + (15)(36)(24) + (16)(35)(24) 
+ (15)(46)(23) + (16)(45)(23) + (25)(36)(14) + (26)(35)(14) 

+ (25)(46)(13) + (26)(45)(13) + (35)(46)(12) + (36)(45)(12)]} 
+(56) [(17)(28)(34) + (18)(27)(34) + (17)(38)(24) + (18)(37)(24) 
+(17)(48)(23) + (18)(47)(23) + (27)(38)(14) + (28)(37)(14) 

+ (27)(48)(13) + (28)(47)(13) + (37)(48)(12) + (38)(47)(12)]} 

+D 422 {(67) [(15)(28)(34) + (18)(25)(34) + (15)(38)(24) + (18)(35)(24) 

+ (15)(48)(23) + (18)(45)(23) + (25)(38)(14) + (28)(35)(14) 

+ (25)(48)(13) + (28)(45)(13) + (35)(48)(12) + (38)(45)(12)]} 
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+ (57) [(16) (28) (34) + (18) (26) (34) + (16) (38) (24) + (18) (36) (24) 

+ (16)(48)(23) + (18)(46)(23) + (26)(38)(14) + (28)(36)(14) 

+ (26)(48)(13) + (28)(46)(13) + (36)(48)(12) + (38)(46)(12)] . . .} 

+^422 {(15)(26)(37)(48) + (15) (26) (38) (47) + (15) (27) (36) (48) 

+(15) (27) (38) (46) + (15) (28) (36) (47) + (15) (28) (37) (46) 

+(16) (25) (37) (48) + (16) (25) (38) (47) + (16) (27) (35) (48) 

+(16) (27) (38) (45) + (16) (28) (35) (47) + (16) (28) (37) (45) . . .} . (B.29) 

If we contract this with (78) we get a six-O integral, which we have 
already calculated. In particular, if we compare the (12) (34) (56) elements 
we find 

NA 422 + 2^422 + 4C 4 22 = ^42, (B.30) 

while comparing the (12) (35) (46) elements gives 

NC 422 + 4^422 + 2^422 = #42- (B.31) 

On the other hand, contracting with (67) gives 5^ = 0. Then looking at 
the coefficient of (12) (34) (58) tells us 

^422 + (N + 1)^422 + 4^422 = 0, (B.32) 

while (12) (35) (48) gives 

2C 42 2 + (N + 2)^422 + 2^422 = 0. (B.33) 

Contracting with (45) also gives zero. The coefficient of (12) (36) (78) 
yields 

^422 + (N + 3)C 42 2 + 2£>422 = 0, (B.34) 

while that of (12) (37) (68) requires 

#422 + ^22 + (N + 4)^422 = 0. (B.35) 

From the last four equations we find 

N 3 + 8iV 2 + 13JV - 2 
2 " N* + 7N + U 54225 



40 



N 2 + 5N + 2 
6422 " W + 7N+U B422 > 
2(N + 2) 

° A22 ~ -Ni + 7N + U Bi22 > 
TV - 2 D 

^ 422 " ~AT2 + 7iV + 14 5422 - 

(B.36) 

Combining these with the preceding equation gives 

N 2 + 7N+U 

422 ~ ~N{N -1){N+ 1){N -2){N + 2){N + 4){N + 6) ' 

(B.37) 

so altogether 

N 3 + 8N 2 + 13N - 2 

A422 '- 



B422 = 
C422 = 
-D422 = 



N(N -1)(N+ 1)(N -2)(N + 2)(N + 4)(7V + 6) ' 

/ N 2 + 7N + 14 \ 

~ \N(N -1)(N + 1)(N -2)(N + 2)(N + 4)(N + 6) J ' 

/ N 2 + 5N + 2 \ 

~ \N(N -1)(N + 1)(N -2)(N + 2)(N + 4)(N + 6) ) ' 
2(N + 2) 



N(N -1)(N+ 1)(N -2)(N + 2)(N + 4) (TV + 6) ' 

N — 2 

rp _ 

422 ~ 7V(7V-l)(7V + l)(7V-2)(7V + 2)(7V + 4)(7V + 6)' 

(B.38) 

Finally, consider the case where there are two sets of four identical indices: 



J fJ,(dR) O aai O aa . 2 O aa3 O aa4 Of3a 5 O '/3a 6 O '/3a 7 0/3(1$ 

= A 44 {(12) (34) (56) (78) + (12) (34) (57) (68) + etc.} 

+B AA {(12)(35)(46)(78) + (12) (35) (47) (68) + etc.} 

+C A4 {(15)(26)(37)(48) + (15) (26) (38) (47) + etc.} (B.39) 



41 



Now contract (78) with the eight-0 integral to find A AA and _B 44 in terms 
of A42 and B42, for a 7^ (3 

(N + 2)A 44 [(12)(34) + (13)(24) + (14)(23)](56) 
+B44 {4[(12)(34) + (13) (24) + (14)(23)](56) 

+(N + 4) [(12) [(36) (45) + (35)(46)] + (13)[(26)(45) + (25)(46)]...]} 
+2C 44 [(12)(35)(46) + (12)(36)(45)...] 
= A 42 [(12)(34)(56) + (13)(24)(56) + (14)(23)(56)] 

+£42 [(12) (35) (46) + (12)(36)(45) + (13)(25)(46) + (13)(26)(45)...], (B.40) 

so 

^42 = (N + 2)A AA + AB AA , 

B A2 = (N + A)B AA + 2C 44 . (B.41) 

We have two equations for our three unknowns, A AA , B AA , C 44 . Two 
additional equations can be obtained by contracting Eq. ( |B.39| ) with (45), 
which gives zero since a 7^ (3: 

[A AA + (N + 5)5 44 ][(12)(36)(78) + (12) (37) (68) + (12) (38) (67) + (13) (26) (78) 
+ (13)(27)(68) + (13)(28)(67) + (16)(23)(78) + (17)(23)(68) + (18)(23)(67)] 

+ [3B AA + (N + 3)C 44 ][(16)(27)(38) + (16) (28) (37) + (17) (26) (38) 

+ (17)(28)(36) + (18)(26)(37) + (18)(27)(36)] = 0. (B.42) 

Thus 

A AA + (N + 5)B AA = 0, 

W AA +(N + 3)C AA = 0. (B.43) 
Combining these with 

B A2 = (N + A)B AA + 2C 44 , (B.44) 

we find 

A (iV + 3)(iV + 5) 

44 (N- 1)N(N + l)(JV + 2)(iV + 4)(JV + 6)' 
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B = -(N + 3) 

44 (N — 1)N(N + 1)(JV + 2)(N + 4)(7V + 6) ' 

c ? 

44 (N- l)N(N + l)(N + 2)(N + A)(N + 6)' 

(B.45) 

C Calculation of Fluctuations 

We explain here how to calculate fluctuations, which require the evaluation 
of 

\M 2 b \ 2 = £ X a X* c X d X: x 

a,b,c,d 

£ J ^dR)\{d5 n WO ab e lS «O aa O lb e-^0 lc 5b e^ 

Qf,/3,7,e « 

(C.l) 

Using the results of Appendix [B], the contraction of X a X*X d X* with the 
integrals over the orthogonal group are found to be 
Category I: 



^ 8 [6(E|A 2 |) 2 +3EA 2 EA* 2 +48|X fc 2 |E|X 2 |+12(X fc 2 EX* 2 +X fc * 2 EA 2 )+24|X 6 2 | 2 ]. 

(C.2) 

Category II: 

A 44 [(£|X 2 |) 2 + 4|X 2 |£|X 2 |+4|X 2 | 2 ] 

+5 44 [3(£|X 2 |) 2 + £X 2 £X* 2 + 36|X 2 |£|X 2 | + 8(X 2 £X* 2 + X* 2 £X 2 ) + 16|X 2 | 2 ] 
+C 44 [2(£|X 2 |) 2 + 2£X 2 £X* 2 + 8|X 2 |£|X 2 | + 4(X 2 £X* 2 + X* 2 £X 2 ) + 4|X 2 | 2 ]. 

(C.3) 

Category III: 



a u [j:x 2 j:x* 2 + 2(X 2 £X* 2 + X* 2 £X 2 ) + 4|X 2 | 2 ] 

+5 44 [2(£|X 2 |) 2 + 2£X 2 £X* 2 + 32|X 2 |E|X 2 | + 10(X 2 £X* 2 + X 6 * 2 £X 2 ) + 16|X 2 | 2 ] 
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+C 44 [4(E|X 2 |) 2 + 16|X 2 |E|X 2 | + 4|X 2 | 2 ]. 

(C.4) 

Category IV: 



A 2222 [|X 2 | 2 ] 

+ J B2 222 [16|X 2 |E|X 2 | + 4(X 2 EX* 2 + x* 2 e^ 2 ) + 8|X 2 | 2 ] 
+C 22 22[4|X 2 |E|X 2 | + (X 2 EX* 2 + X* 2 E^ 2 ) + 6|X 2 | 2 ] 

+ j d 2222 [2(e|x 2 |) 2 + 4|x 2 |E|x 2 | + (x 2 ex* 2 + xi 2 Y,x 2 ) + 3|x 2 | 2 + E^ 2 E^* 2 ] 

+E 2222 [4(E|X 2 |) 2 + 24|X 2 |E|X 2 | + 6(X 2 EX* 2 + x* 2 E^ 2 ) + 6|X 2 | 2 + 2EX 2 EX* 2 ]. 

(C.5) 

Category V: 



A 422 [2|X 2 | 2 + |X 2 |E|X 2 |] 
+5 422 [3|X 2 |E|X 2 | + (E|x 2 |) 2 + 2|X 2 | 2 ] 
+c 422 [io|x 2 |E|x 2 | + 3(x 2 ex* 2 + x; 2 j:x 2 ) + 8|X 2 | 2 ] 

+£422[3(E|^ 2 |) 2 + 26|X 2 |E|X 2 | + 5(X 2 EX* 2 + x* 2 E^ 2 ) + 8|X 2 | 2 + E^ 2 E^* 2 ] 
+^42 2 [2(E|X 2 |) 2 + 8|X 2 |E|X 2 | + 4(X 2 EX* 2 + xi 2 Y,x 2 ) + 4|X 2 | 2 + 2EX 2 EX* 2 ]. 

(C.6) 

Category VI: 

Half the terms have 

A 422 [2|X 2 | 2 + X 2 EX* 2 ] 

+5 422 [EV 2 ex* 2 + 2X* 2 e^ 2 + x 2 ex* 2 + 2\x 2 \ 2 ] 

+C 422 [12|X 2 |E|X 2 | + 4X 2 EX* 2 + 8|X 2 | 2 ] 

+A22[2(E|X 2 |) 2 + 20|X 2 |E|X 2 | + 6X 2 EX* 2 + 10X* 2 E^ 2 + 8|X 2 | 2 + 2EX 2 EX* 2 ] 
+i?422[4(E|X 2 |) 2 + 16|X 2 |E|X 2 | + 4|X 2 | 2 ], 

(C.7) 

while the other half have the complex conjugate result. 
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Category VII: 
Half the terms have 



A 62 [(3X; 2 T,X 2 + 6|X 2 |£|X 2 | + 6|X 2 | 2 ] 

+B 62 [6{j:\X 2 \) 2 + 3£X 2 £X* 2 + 42|X 2 |£|X 2 | + 12X 2 £X* 2 + 9X b * 2 £X 2 + 18|X, 



'2 1 21 



(C. 



while the other half have the complex conjugate result. 

We can use these results to determine the fluctuations and the direct 
CP violation, Eq. QS.3| ). In the latter instance, we see that the factor 
e r(5 a -(5 7 )^ e 2(5 a -5 6 ) _ e -t(s s -6 e )^ interchanges the second line of Category II in 

Table ( |5TT| ) with Category III, and the second and third lines of Category V 
with the two lines of Category VI. The result is 



<(\Ml\-\M;\) > 



2 (A 



11 



B 



li 



2C 



i-i ) 



X 4 |X 2 | -X 4 X 2 X* 2 



+4x 2 |x b 2 ||x 2 | - 2N 2 (x b 2 x* 2 + x; 2 x 2 ) 



'1- < cos 25 >' 



+ [N 2 (A 422 + 3B 422 - 2C 422 + 6D 422 - 8E 422 ) [2 {X^JX 2 ] - (X 2 X^ + XfX 2 )} 

+N\B 422 + D 422 - 2E 422 )(]X 2 \-X 2 X^ 2 )] < cos<5 > 2 (1- < cos 26 >) 

(C.9) 

In the limit of large N, the most important coefficients are A 44 and A 422 , 
each of which is approximately N~ A . Thus for large N, 



< \Ml\-\Ml\) > 



2 j(jX 2 | - X 2 X* 2 + A( 2 |X 2 ||X 2 | - X 2 X* 2 - X* 2 X 2 )) (1- < cos 25 > 2 ) 



+ — (2|X 2 ||X 2 | - X 2 X* 2 - X* 2 X 2 ) < cos 5 > 2 (1- < cos 25 >] 



(CIO) 
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D Correlations in the Dyson Model 

We further define n- level cluster functions, T n , via 

oo . „ 27r 

£/•••/ Ru(S 1 ,...8 n )A(8 1 )---A(8 n )d8 1 ---d8 n 

n=0 J J ° 

00 ( — l)™- 1 r r2n 

= exp]T l , / ••• / T n (8 1 ,...8 n )A(8 1 )---A(8 n )d8 1 ---d8 n , 

(D.l) 

where A (8) is an arbitrary function. In particular, 

R2{8i,8 2 )=T l {8 1 )T l {8 2 )-T 2 {8 1 ,8 2 ). (D.2) 

In a tour de force, Dyson calculated T 2 and showed how to calculate, 
in principle, the higher-order Ts. Explicitly, for even N = 2m, and with 
p,q = -m + \,...m—\ 



where e(C) = C/ICI- 
It is clear that 



/ 2 " [ 2 * d8d<f)R 2 (8,(f)) = N(N - 1), 

JO JO 

/ d5Ti(5) = AT, 
Jo 

/*27T /"27T 

/ / d8d(pT 2 (8, 0) = AT. 
jo jo 



We need to calculate 

r 2n r 2n 



(D.4) 



/ W / * dSd<f)T 2 {S, 0) cos {8 - 0). (D.5) 

JO JO 
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We first compute, for a half-integer v 



/ / tW#e(<5 - (fieW-V = — , (D.6) 
jo jo f 



from which we ascertain 



p p \ 1 1 



m 



= -3EH L r + Ji T = ~ 2 ™ + ^ + 



For 7^ an integer, we find 

r-2-TT /-27T 



/ / d9d<j)e l ^ e ~^ = 0, (D.8) 
jo jo 

so the only contributions to the second part of our calculation come from 
p — q ± 1 = and 

dfcty£ ( 2 + - + -1 e^ 9 -^ cos(d - 0) 
8tH jo jo ^ \ q Pj 

= - V 2 + — ^- + F - + 2 + + F - 

4^ LV p+1 p p-l p 

U, 1 1 . 

l r 1 1 1 1 , 

= 2-2m+-[ r r + r + 



4 m + \ m — ^ —m — \ —m + \ 

= 27V- [—*— + — *—]. (D.9) 
l N+l N-V y ' 

Altogether, then, 

r2ir p2ir 1 1 

^ J Q ^0T 2 (^,0)cos(^-0)=iV + FTT + ^- T . (D.10) 
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It follows that 

[l[d8 k V = 1 

J k 

= _^i- Iy[ i V+F L T + ^L T ]. ( D.ll) 



k 

Assembling this with our previous calculation gives 



< \M b \ 2 >= ' 



N 2 + 1 



n + 2 L 2 ^ 2 + - (iv.-iKjy-D ^W - Ar|x2|) 

(D.12) 

Let us now approximate the correlation functions described by Dyson. 
The correlations are of a very short distance, of order 0(1/ N). Thus it 
makes sense to use the approximation 

T 2 (9 1 ,9 2 ) = ^6(9 1 -9 2 ), (D.13) 
where the normalization is determined by the relation 

/ d9T 1 (9) = N. (D.14) 
Jo 

We proceed analogously writing 

t-2-K rl-K rl-K 



rZTT rZTT pZTT 

/ dOj d9 2 d9 3 R 3 (9 1 ,9 2 ,9 3 ) = N(N-l)(N-2), 
Jo Jo Jo 

/■27T r2n e-l'K 

j d9 1 J Q d9 2 J Q d9 3 [T 3 (9 1: 9 2: 9 3 )-T 1 (9 1 )T 2 (9 2: 9 3 )-T 1 (9 2 )T 2 (9 u 9 3 ) 
-T^T^e*) +T 1 {9 l )T l {9 2 )T l {9 3 )\ , (D.15) 
from which we conclude 



T 3 (9 ± , 9 2 , 9 3 ) = -5(9, - 9 2 )5(9 2 - 9 3 ). (D.16) 
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Similarly 



Ra(9i, 9 2 , 9 3 , 4 ) = N(N - 1)(N - 2)(N - 3) 

x [-T 4 (9 U 9 2 , 9 3: 9 A ) +T 3 T 1 ... + T 2 T 2 ... - T 2 T{T X ... + , 

(D.17) 
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where the arguments are implicit. In this way we determine the approxima- 
tion 



74(0i, 2 , 3 , 4 ) = —5(6 1 - 9 2 )5{9 2 - 9 3 )5(9 3 - 4 ). 

7T 

The required integrations over the eigenphases are 



(D.18) 



o Jo 2 N(N-l) 

2w fl-K 



JO JO 



/•27T /"27T 

« / / d9 x d9 2 -^ 
Jo Jo 



N(N - 1) 



N(N-1) 



N-l 



(D.19) 



and 



and 



27T /-27T /■27T 



j fl j fl j fl ^3(01,02,03) f ,t(20 1 -0 2 -g 3 ) 

n 70 ./n dBl<W2d03 N(N-l)(N-2) e 

2n r 1-K fl-K 



> ' ' ' <I0\<W 2 <W A 
io Jo Jo 



x 



at p i(2ei-e 2 -e 3 ) 
2iV 



N(N-l)(N-2)' 



(D.20) 



/o Jo Jo Jo 123 4 jV(iV- l)(7V-2)(7V-3) 

/•27T /-27T /■27T /-27T 



-02+0 3 -0 4 ) 



x 



/o jO jo jo 
3iV 



d9 x d9 2 d9 3 d9 A 



{ 5(0i - 2 )5(0 2 - 3 )5(0 3 - 4 ) 

I 7T 
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+ (^) 2 [<K#i - W(# 3 - 04) + 8(0i ~ 3 )S(9 2 - 4 ) + - #4)5(# 2 - # 3 )} 



3 i(fl-02+03-f4) 



X 



X(X-l)(X-2)(X-3) 
6X + 2X 2 



X(X- l)(X-2)(X-3)' 



(D.2i; 



Of course, we have approximated the correlations and should only work to 
leading order in X. We note that in the various tensors, the leading piece is, 
in every case, the A term, which is, to leading order, 1/X 4 . We collect the 



results in Table D.l 





Tensor 


Terms 


Eigenphase 
Integral 


I 


6(£|X 2 |) 2 + 3£X 2 £X* 2 + 48|X 2 |£|X 2 | 
+12(X b 2 £X* 2 + X 6 * 2 £X 2 ) + 24|X 6 2 | 2 


X 


1 


II 


(E|* 2 |) 2 + 4|X 2 |£|X 2 |+4|X 2 | 2 


2X 2 


1 


III 


EX 2 £X* 2 + 2(X 6 2 EX* 2 + XfJ2X 2 ) + 4|X fe 2 | 2 


X 2 


-1/X 


IV 


l* 6 T 


X 4 


2/X 2 


V 


2\X 2 \ 2 + \X 2 \E\X 2 \ 


4X 3 


-1/X 


VI 


4|X fe 2 | 2 + X^X* 2 + XfY.X 2 


X 3 





VII 


3Xf EX 2 + 6|X 2 |E|X 2 | + 6|X 2 | 2 


2X(X- 1) 
2X(X-1) 


-1/X 



Table D.l: Leading terms in X for the fluctuations in the amplitude squared 
in Dyson's model for the various categories. 

Combining these results and taking only the pieces leading X we find 



<|A^I 2 >=^(E|X 2 | + |X 6 | 



(D.22) 
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